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Abstract 


The Sherman-Morrison-Woodbury formula-based algorithm is used in this 
work to enhance the computational performance of the characteristic basis 
function method (CBFM) when applied in a context of 3D full-wave model to 
the scattering by complex-geometry precipitation particles. The improvement 
brought by this algorithm, particularly to the compression rate achieved by the 
CBFM, enables us to simulate electrically larger particles while maintaining a 
Satisfactory level of accuracy and a reasonable computational cost. 





Index Terms: Scattering by hydrometeors, Volume Integral Equation Method 
(VIEM), Method of Moments (MoM), Characteristic Basis Function Method 
(CBFM), Sherman-Morrison-Woodbury formula-based algorithm, Adaptive 
Cross Approximation (ACA), Discrete Dipole Approximation (DDA), DDScat. 
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e VIEM/MoM-based 3D full wave model for scattering by complex- 
shaped precipitation particles. 


e Application of the domain decomposition-based CBFM 


2) Efficient generation of the CBFs 


e Impact of the size of the blocks on the compression rate 
e Sparse representation of the MBFs 
e Sherman Morrison Woodbury Formula-based algorithm 





3) Numerical Analysis 





e Numerical Analysis : NESCoP vs DDSCat 


e Impact of the efficient calculation of the CBFs on the overall 
computational performance of NESCoP 
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Scattering by large and complex-shaped 
precipitation particles 





Calculate the single-scattering properties of 
particles with complex geometries modeling 
snowflakes of various sizes and sha 











Pristine crystals (a) simulated using the Snowflake 
algorithm and aggregate snow particles (b) 














Scattering by large and complex-shaped 
precipitation particles 









Calculate the single-scattering properties of 
particles with complex geometries modeling 
snowflakes of various sizes and shapes. 














Discrete dipole approximation (DDA) 
representation of a sphere [1] 


Discrete Dipole approximation 
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= m : complex refractive index 
= k: wavelength number 
= d: grid spacing 


11] B.T. Draine and P.J. Flatau, “Discrete-dipole approximation for 
scattering calculations,” in JOSA A, 11(4), 1994, pp. 1491-1499. 














Discrete dipole approximation (DDA) 
representation of a sphere [1] 








Scattering by large and complex-shaped 
precipitation particles 











Calculate the single-scattering properties of 
particles with complex geometries modeling 
snowflakes of various sizes and sha 


Discrete Dipole approximation 
(DDA) : DDScat, ADDA, ... 
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Orientation averaging 


11] B.T. Draine and P.J. Flatau, “Discrete-dipole approximation for 
scattering calculations,” in JOSA A, 11(4), 1994, pp. 1491-1499. 










Scattering by large and complex-shaped 
precipitation particles 






Calculate the single-scattering properties of 
particles with complex geometries modeling 
snowflakes of various sizes and sha 











Discrete Dipole approximation 
(DDA) : DDScat, ADDA, ... 


Characteristic Basis Function Method 
(Direct Solver-based method) 
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VIEM/MoM-based 3D full-wave model 


3D full-wave model, based on the volume integral equation method (VIEM) 








= Frequencies of interest : 
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VIEM/MoM-based 3D full-wave model 


3D full-wave model, based on the volume integral equation method (VIEM) 
Frequencies of interest : 
(13 - 200 GHz) 


Pristine crystal particle 
of effective radius 
a, = 0.714476318 mm 
discretized into Nb, = 
12222 elementary cubic 
cells 








Method of Moments 


The particle is discretized into N 
cubic cells Q, ,of side c, 
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Z is the 3Nb, x 3Nb, full 
matrix representing the 
interactions between the 
different cells. 





Use of piecewise 
constant basis 
functions 








VIEM/MoM-based 3D full-wave model 


3D full-wave model, based on the volume integral equation method (VIEM) 
Frequencies of interest : 
(13 - 200 GHz) 


Pristine crystal particle 
of effective radius 
a, = 0.714476318 mm 
discretized into Nb, = 
12222 elementary cubic 
cells 






Method of Moments 


The particle is discretized into N 
cubic cells Q, ,of side c, 


< ls 2 
ie === 
^ .10' " JRe(e,) 


Z is the 3Nb, x 3Nb, full 
matrix representing the 
interactions between the 
different cells. 












Example: M=4 
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Example : M - 4 
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Application of the domain decomposition-based CBFM 


3D full-wave model, based on the volume integral equation method (VIEM) 
Frequencies of interest : 
(13 - 200 GHz) 


Pristine crystal particle 
of effective radius 
a, = 0.714476318 mm 
discretized into Nb, = 
12222 elementary cubic 
cells 










Method of Moments 


The particle is discretized into N 
cubic cells Q, ,of side c, 
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Impact of the size of the blocks on the compression rate 





Effect of the size of the CBFM blocks on the Compression rate 
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100 : 200 
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Nb & Nbe 
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x Number of basis functions per block T id 
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Enhancement Techniques : Sparse Representation of the MBFs 
VAL E! 55 FE! 
MBFs~ “IPWs 
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= Generation of the CBFs 


fs, 1s a threshold factor used to down-select 
the elements of Z whose magnitudes are 
significant compared to | ZÎ , ,l. 
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Sparse Direct Solver | 
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X. Chen, & al, "Accelerated Direct Solution of EM 
Scattering via CBFM With Sherman-Morrison-Woodbury 
Formula-Based Algorithm," IEEE Transactions on Antennas 
and Propagation, Oct. 2016. 





Application of a multilevel Sherman- 
Morisson-Woodburry Formula- 
based algorithm (SMWA) 


The Sherman-Morrison Formula provides 
an explicit formula of the inverse of 
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Application of a multilevel Sherman- 
Morisson-Woodburry Formula- 
based algorithm (SMWA) 





1-Level SMWA : 
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X. Chen, & al, "Accelerated Direct Solution of EM 1 0 + 0 
Scattering via CBFM With Sherman-Morrison-Woodbury 0 1 U 0 
Formula-Based Algorithm," IEEE Transactions on Antennas 21 


and Propagation, Oct. 2016. 
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Application of a multilevel Sherman- 
Morisson-Woodburry Formula- 
based algorithm (SMWA) 





1-Level SMWA : 
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X. Chen, & al, "Accelerated Direct Solution of EM 
Scattering via CBFM With Sherman-Morrison-Woodbury 
Formula-Based Algorithm," IEEE Transactions on Antennas 
and Propagation, Oct. 2016. 








Enhancement Techniques : Sherman-Morrison-Woodbury Formula 
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3-Levels SMWA : 







Application of a 
multilevel Sherman- 
Morisson-Woodburry 
Formula-based 
algorithm (SMWA) to 
solve : 
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Zhang, Yang, et al. "Fast Direct Solution of 
Composite Conducting-Dielectric Arrays 
Using Sherman-Morrison-Woodbury 
Algorithm." Progress In Electromagnetics 
Research M 49 (2016) 
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Numerical Analysis : NESCoP vs DDScat 


N ES CoP Numerically Efficient Scattering by Complex Particles 
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Numerical Analysis : NESCoP vs DDScat 


Orientational averaging 
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Numerical Analysis : NESCoP vs DDScat 


Impact of the 
SMWF-based 


algorithm on the 
computational cost 
of NESCoP 








@ f = 150 GHz over 2700 id Snow aggregate a0013 1s divided into 27 and 4 blocks 
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Numerical Analysis : NESCoP vs DDScat 


Impact of the 
SMWF-based 


algorithm on the 
accuracy 
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Numerical Analysis : NESCoP vs DDScat 


a, = 1.61 mm; 
d,,=11.45 mm; 
16 xf x 200 GHz 
0.17 € x, > 4 
Imlkd < 0.37 
Nb, =140896 cells 


Impact of the SR 
of the MBFs on the 


computational cost 
of NESCoP 








TABLE III: Relative difference E,.;(%) and computational proformance (CPU time in minutes, size of Z° and compression 
rate) when calculating averaged scattering by a0006 with 2701 id at f = 150 GHz, depending on the use or not of the DR 
and SR techniques, on hp and fsr. < Nb, > refers to the average extended block size and p; sr (%) is the stored proportion 
(in percentage) of the sparse self-coupling MoM sub-matrix Z^ 
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Numerical Analysis : NESCoP vs DDScat 
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Numerical Analysis : NESCoP vs DDScat 
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TABLE : CPU time (min) needed by NESCOP, with 1L and 
2L CBFM-E, and DDScat to calculate averaged scattering by 
a0006 with 703 and 2701 id and 900 and 1800 to for 16.67 
GHz < f < 200 GHz. 


۲ 10/10 | 703 id | 900 to | 1800 to | 2701 1d 
IL CBFM-E | 126.72 | 527.67 | - | - [| 1706 


JL CBEM-E | 32947 | 53215 | — | — [10548 
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Conclusions and Perspectives 





3D full wave model comparable to the DDA 
in terms of accuracy when providing higher 


computational performance particularly 
with orientational averaging of the EM 
scattering due to a high CR achieved by the 
CBFM. 


* Exploring further techniques for 
efficient calculation of the CBFs 

** MPI parallelization of the codes 

* Implementation of the wideband 
code for the CBFM 
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